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We formulate lattice fermions in a way that encompasses Wilson fermions as well as the static and non- 
relativistic approximations. In particular, we treat m^a systematically (m, is the fermion mass) showing how to 
understand the Wilson action as an effective action for systems with p <C ruq. The results show how to extract 
matrix elements and the spectrum from simulations with mqa ~ 1, which is relevant for the charm quark. 



1. INTRODUCTION 

Current numerical calculations of lattice QCD 
have lattice spacings such that niqa <^ 1 for 
charm and bottom quarks. To extract meaning- 
ful results one must understand how to formu- 
late fermions under these circumstances. One 
approach is to formulate this problem starting 
from the large mass limit. Our approach is to de- 
duce an action that contains the Wilson ^ , static 
m , and non-relativistic [|| formulations as special 
cases. An important question is how much of the 
flexibility of our formulation should be exploited 
to study, for example, quarkonia and heavy-light 
systems. A related and practical issue is how one 
can interpret numerical results obtained with the 
Wilson action when m^a « 1. This paper offers 
a brief discussion of these issues. 

The key idea is that one should analyze every 
heavy-quark theory, including the Wilson theory, 
as an effective (cutoff) theory. After all, every 
numerical simulation has a 7^ 0. In other words, 
one should avoid being constrained by an implicit 
desire to take the continuum limit by brute force. 
Naturally, the goal remains the same, to extract 
universal (i.e. cutoff-free) information, but one 
would like to do so as efficiently as possible. 

With any effective theory, one must understand 
its limitations. It is worth recalling that the 
characteristic momenta of the systems that we 
would like to study are much smaller than niq. 
In quarkonia the typical momentum scales are 
as{mq)mq and as{mq)^mq, so pa !C 1 for a wide 



cal momentum is roughly Aqcd, independent of 
uiq. Below we argue that the Wilson action de- 
scribes heavy quarkonia with the same validity as 
the non-relativistic effective theory and describes 
heavy-light systems with the same validity as the 
non-relativistic and static effective theories. 



2. THE NEW ACTION 

The familiar heavy-quark effective theories 
treat time and space asymmetrically. Therefore, 
it should not be too surprising that our result jj] 
boils down to introducing a second hopping pa- 
rameter to the familiar Wilson action. 
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In heavy-light systems the typi- 
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This action describes massless quarks when kq = 
Kg ^ 1/8. It also describes static quarks when 
for Ks = 0. Intuitively, there ought to be a trajec- 
tory in the kq-ks plane that describes somewhat 
massive quarks. 

Eq. (|l]) actually describes the simplest in a 
class of actions, in which lattice artifacts are cor- 
rected by the Symanzik improvement program. 
Further details appear elsewhere [||, but setting 
'^0 = '^s = '^j the most important term is the 
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"clover" term 



(2) 



where C^i, is a lattice approximant to F^^. This is 
the form usually used in numerical work, but for 
perturbative calculations and contact with con- 
tinuum theories it is convenient to write n = x/a, 
(2k)"1 = 4 + amo, and (2^)^/2^^ ^ o?/'^^{na). 

3. THE CONTINUUM 

This section is in Minkowski space; everything 
else is in Euclidean space. The continuum ac- 
tion is 5 = / d'^x'ilj{ip + m)tp. Choose the gauge 
Aq = 0, split the fermion into two-component 
fields, and redefine the fields (p i-> e™''*(/) and 
X 1-^ e~*'"°*x. Then, using the appendix, 
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+(j)W-(TX^ -XD-(T(I)) . 



Because x appears quadratically in the functional 
integral, it can be integrated out, leaving an ef- 
fective theory for (j) with action 

f d^xcj)^ (ido-D-cT ] D-ct]cI> (4) 
J \ 2m + ido ) 

For systems for which idoD-crcf) <C mD-crcf) the 
non-local expression can be expanded 
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Higher time derivatives can be eliminated using 
the equations of motion (e.g. [9o,-D] = E) ||], 
(D-cr)2 = +ia B, etc. 

Usually only the D'^/{2m) and the a-B/{2m) 
terms are treated explicitly, with short-distance 
effects of the higher terms compensated by ad- 
justing the couplings of these two terms. Overall, 
this type of analysis is very successful, for exam- 
ple, motivating static potential models and devel- 
oping of heavy-quark symmetry. One limitation 
is obvious: the field (p now creates states with 
energy-momentum relation E = p^/(2m) instead 
oi E = m +p^/(2m). Hence the splitting be- 
tween sectors with a different number of quarks 



does not come out right. Usually they can be 
corrected by hand, but one should keep an open 
mind to more subtle dynamical effects that might 
alter, say, mixing between heavy-quark mesons 
and glueballs. 

4. THE LATTICE 

Suppose that a lattice theory has the following 
1/M expansion for the Hamiltonian: 



H = a) 



Ml - 



2M2 



ia-B 
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-|- anti-quark terms, 

but the M's are not all the same. (Naturally, in a 
lattice theory the Hamiltonian is defined through 
the transfer matrix formalism.) They play the 
following roles: 

• Ml counts quarks and anti-quarks * 

• M2 fixes spin- averaged sphttings ***** 

• M3 fixes hyperfine structure *** 

• M4 fixes fine structure ** 

Of these "masses" M2 is the most important 
physically. As in sect. ^ (where Mi — and 
M2 = m) the fact that Mi ^ M2 can be cor- 
rected by hand. If that strategy is impractical or 
distasteful, one simply switches back to eq. (|^); 
sect. shows how the k's can be tuned so that 
Ml = M2. Similarly, the other masses multiply 
operators that are easy to interpret physically, so 
their deviations could be taken into account as 
perturbations. Alternatively, they can be cor- 
rected as part of the Symanzik improvement of 
eq. i§. 

All actions in the class exemplified by eq. (|l|) 
take the form of eq. (^, in the limit \p\ ^ M2, 
pa <C 1 , even without restriction on aM2 . For 
example, at tree level the Wilson action has 

Ml ^ a-Hog{l + moa), (7) 

2 + 4moa + m^a^ 
2mo{l + moa){2 + moa) (g) 

(2mo)-i, 



(2M2)- 
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and 

(2M3)-i = (2M2)-i - 



1 - c 



2too(1 + moa) 



(9) 



These expressions are corrected at 0(.go). The 
first two are easily extracted from the (free) 
energy-momentum relation 



cosn(iiaj = 1 + 77 



1.2^2,2 ^^2^2 



2 1 + moa + ia'^p^ 
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Note that (2Af2)"^ = dE/dp^ arises from p'^ and 
the Wilson term, whereas the chromomagnetic 
term arises only from If)^ . This leads to a mis- 
match with (2M3)~^, which is removed by the 
clover term . 

5. A TOY MODEL 

To derive eq. (H) one must set up the transfer 
matrix formalism. First one separates the four- 
component fields into two-component fields, using 
the appendix. The manipulations become quite 
cumbersome [Q, so consider instead a toy model 
with the most important properties: 



S = ^(t)\{d^ + m)(t)t + x\{do + m)xt 
t 

+ip4>\x\ - 'iPXt<l>t- 
Then VipVip = Ht d(t>ld(t>tdxldxt and 

e'S ^ J^g~(l+m)(0l0t+xlxt)x 

t 

YlT{4+i,xl+i;(t>t,Xt), 
t 

where 



(11) 



(12) 



T{(f)\x'']<i), x) = e-'J"^ X e"^ '^+>^^e'P>^'^. (13) 

The rules of Grassman integration imply 
{a, a)} = 1 = {b,b^}, where 



(14) 



b = \/l + mx 



are operators in Hilbert space [Q. The Fock vac- 
uum satisfies a|0) = b\Q) 0. The other states 
are \q) = d^O), \q) = b^O), and \qq) = d^b^O). 



Up to Jacobian factors, the matrix elements of 
the transfer matrix are the coefficients of mono- 
mials in (Grassman numbers) a, a\ 6, and 6^, 
when T{(p\x^^ 4'i x) is expressed as a polynomial. 
Thus {q\f\q) = {q\f\q) = 1 -l- m, and in the neu- 
tral sector 



mi) 



-ip{l + m) l+p^ 



(15) 



and all other elements are zero. This model can 
be solved exactly, by diagonalizing (z|T|j). Upon 
reconstructing the operator form and expanding 
in p, one obtains the analog of eq. (|^) . 

6. THE NORMALIZATION 

and b^ create normalized states, 0^ and 
do not. Applying eq. (nih to Wilson fermions, 
with m amo — ^a^A^^one sees that 

/ \ 1/2 

a(a;) = n + amo - ifl^A^^' j tpix) (16) 
is canonically normalized; in momentum space. 



(17) 



1 /2 

a{p) = (1 + amo + ia^) V'(p)- 

At p = this factor is e*^i/^. 

This fact has been noticed in several phe- 
nomenological analyses |^,||, which suggested 
using e'^/'^Q{x)^^j,q{x) and similar bilinears to 
determine heavy-light decay constants, using 
Ansatze such as e*^ = 1 + 1/(2k) - 1/(2kc). Com- 
bining the transfer-matrix analysis with mean 
field theory suggests a different Ansatz. Including 
gauge fields and using the notation of eq. 

a„ = (1 - «;B)y^i/-™, (18) 

has unit normalization ||^, where -BnmV'm = 
Y.i Un,itpn+i + Ul_.^^il;n-i- In mean-field theory 
one replaces the dynamical C/'s by their average. 
Hence (1 — GkY^'^ipni where k = fi(u) ~ k/Sk^ 
||lC| ], is expected to be correctly normalized to 
high accuracy. As shown in fig. ^ the mean-field 
Ansatz agrees perfectly with numerical data for 
the local current. 

7. MORE ON THE NEW ACTION 

This section concludes with two short com- 
ments on eq. (pi). It is convenient to define 
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auiQ = (1 — 2ko — 6ks)/(2ko) and <^ — ks/kq. 
Then eq. (0) and the denominator of eq. are 
the same, but the numerator of eq. (0) becomes 



A. FOUR TO TWO COMPONENTS 



C(2C + 2moa(l + C) +TO^a^). 



(19) 



An appropriate choice of C imposes Mi = M2. 

The asymmetric action is also natural from the 
point of view of renormalization. Imagine start- 
ing with the Wilson action on a fine lattice, and 
imagine performing a "block spin" to a coarse lat- 
tice. In perturbation theory this entails calculat- 
ing the self-energy S(p). For a massive quark, 
the most natural subtraction point is on shell, 
e.g. E = Ml, p = 0. On the lattice, however. 
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on shell 



^ dp. 



(20) 



on shell 



In the hopping-parameter notation this difference 
implies kq 7^ Kg for the effective action on the 
coarse lattice. 
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Figure 1. {J/ip\V4\J/ip), where V4 = '07o'0 is the 
(unnormalized) local charge. The lower line is the 
mean-field Ansatz Zy/ (1 — 6k) and the upper line 
is Zv/{2k). (The latter form is the one suggested 
by the massless limit.) To get the charge (should 
be 1) divide the data by the curve. 



The Dirac matrices are taken to be 
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In Minkowski space 77 = —i; in Euclidean space 
77 = 1. Fermion fields are decomposed into two- 
component form as follows: 



^ = (0t 



x). 



(22) 
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